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Abstract 

We investigate the aoi^^-ftr~-vertex in the framework of hght cone QCD 
sum rules. We estimate the coupling constant gaoK+K- which is an essential 
ingredient in the analysis of physical processes involving ao(980) meson. Our 
result is somewhat larger than the previous determinations of this coupling 
constant. 



PACS numbers: 12.38.Lg;13.25.Jx;14.40.Cs 



Typeset using REVT^ 



*agokalp@metu.edu.tr 
^oyilmaz@metu. edu.tr 



The structure of light scalar mesons /o(980) and ao(980) has been a controversial problem 
in hadron spectroscopy. In the naive quark model qq [1], ao(980) can be interpreted as 
ao = {uu — dd)/\/2. On the other hand, the strong coupling of /o(980) to kaons suggests 
the structure /o = ss for this state. However, then the almost exact degeneracy of the 
masses of ao(980) and /o(980) cannot be explained. In order to explain the properties of 
these states several proposals have been put forward over the years. A four quark q^lf state 
interpretation with symbolic quark structure /q = ss{v/u, + dd)/^/2 and (Iq = ss{vaI — dd) / \/2 
was proposed in the framework of MIT-bag model where the scalar meson states are spatially 
compact [2]. Another possibility about the structure of ao(980) and /o(980) was suggested 
where these meson states are considered to be bound states of hadrons. This possibihty 
is referred to as their being KK molecules in which case they are considered as spatially 
extended objects [3]. Furthermore, some analyses suggest the qualitative picture that these 
scalar meson states have a compact q^q^ structure that spends some part of its lifetime in 
the KK meson system [4]. 

The radiative decays of 0(1020) meson — > qqJ and (f) — > /o7 provide important tests 
to distinguish among the different possibilities about the structure of ao(980) and /o(980) 
scalar meson states [4,5]. It is generally agreed that the experimental data supports the 
kaon loop mechanism for these decays in both the g^g^ state and KK molecule models 
where these radiative decays proceed by photon emission from an intermediate K^K~ loop 
[5]. Moreover, the study of the reactions (f) — > 7r°7r°7 and — > t^^t]^ and of the interference 
patterns in these reactions have been used to develop arguments about the structure of 
ao(980) and /o(980) states [6]. In the analyses involving ao(980) meson the strong coupling 
constant gaoK+K- plays an important role. 

In this work, we estimate the coupling constant gaoK+K- by employing light cone QCD 
sum rules. This method has been applied to study hadronic properties and in particular it 
has been used for the calculation of hadronic coupling constants [7] . 

In order to study the a^K^ K~ -veitex and to estimate the coupling constant gaoK+K- 
we consider the two-point correlation function 

q)=^j d'xe^^-' < X+(g)|T{j^(x)r(0)}|0 > (1) 

with p and the four-momentum and the interpolating current for K~ meson, the 
interpolating current for meson, and q the four-momentum of K^ state. The interpolating 
quark currents are the axial vector jj^ = u^JuJ^s and the scalar j"" = (uu — dd)/2 current. 
The scalar current is assumed to have a non-vanishing matrix element between the 
vacuum and ao(980) meson state, < oq | | >= Aag, where is called the overlap 
amplitude which was determined by QCD sum rules method [8], and this particular choice 
for the current in terms of quark fields does not imply the pure {uu — dd) / \f2 structure for 
ao(980) meson. 

The correlation function can be written in terms of two independent invariant functions 
Ti and T2 as 

r^(p, q) = iT, ((p + qf, p^) p^ + T2 ((p + g)^ p') q^ . (2) 

We consider the invariant function Ti. In order to construct the theoretical part of the sum 
rule for the coupling constant gaoK+K- we calculate the function Ti in terms of QCD degrees 
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of freedom by evaluating the correlation function in the deep Euclidean region where and 
{p + q)'^ are large and negative as an expansion near the light cone = 0. This expansion 
involves matrix elements of non-local composite operators between kaon and vacuum states 
which defines kaon distribution amplitudes of increasing twist. We retain terms up to twist 
four accuracy since higher twist amplitudes are known to make a small contribution [9]. In 
our calculation we use the full light propagator with both perturbative and nonperturbative 
contributions which is given as [10] 



iS{x, 0) = < 0\T{q{x)q{0)}\0 > 

.2 - ^ 



/t < qq > 



-^^^16^ r rf^ - Uu^^G^^{ux)^)^ + ... . (3) 

After a straightforward computation and performing the Fourier transforms we obtain 

f^Mj, f\ , 2\p^{\-u)q\-q 



Tlip , [p + q) ) = 7-77 r / du^^K{u) — — - 

+ 2-^ Jo 



/ dv / Da,^sK{a.)——- ^ -—^(2^ - 1) . (4) 



{[p+ (1 -tti -vas)q]'^y 



In this expression the functions (p^^K and (ppx are the twist 3 quark-antiquark kaon distri- 
bution amplitudes defined by the matrix elements [11] 

< K{q)\u{x)i^5sm^ >= fKl^K C due'^^'^-^ippKiu) , (5) 

Jo 

< K{q)\u{x)a,,^,siO)\0 >= (l - ^] {q,x, - x,q,) /' due^'^'^-^^^K^u) , (6) 

where — Mj^/rris is the twist 3 distribution amplitude normalization factor and we put 
TTT'u — rrid — 0. We work in the gauge x^A^^ — 0, consequently the path-ordered gauge factor 
is not included in the matrix elements. The twist 3 quark-antiquark-gluon kaon distribution 
amplitude ip^K is defined as [11] 

< K{q)\u{x)aafsgsGiiu{vx)s{0)\0 > = ifax [{qnqaQ'^li - Qi^QaQnf}) - {q^lqpgua - quqiBQua)] 

X j L>Q;,(^3;f(Q;,)e^"^-^("i+""^) (7) 

where Dai = daida2da36{l — ai — 0:2 — cts)- After performing the Borel transformation with 
respect to the variables = —{p + qY and Ql = —p^, we obtain the theoretical expression 
for the invariant function in the form 

T,{Ml Ml) = ^^J^ i-^M^o) + g<^;^M] 

H / dai ip3K{ai,l- Oil- as,a3) [2 1 (8) 

z Jo Juo-ai as \ as / 
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where and are the Borel parameters and 



We hke to note that if we multiply the correlation function by the four-momentum p, we 
obtain 



+6{xo)<K+{q)\\j^{x), rm\0> , (9) 



where the second term results from the differentiation of the function 9{xo) in the T-product 
of the currents. In the 811(3)/^ limit d'^j^lx) = 0, thus the first term on the right hand side 
of Eq. (9) vanishes. The second term can be calculated using the standard commutation 
relations, yielding for the correlation function the Ward identity 

Similar Ward identities were considered in [11] where they were used to obtain relations 
between various pion distribution amplitudes. 

Two-point correlation function satisfies a dispersion relation, therefore we can represent 
the invariant function as 

T.{i,^.f./)-ll isis' . (10) 

We saturate this dispersion relation by inserting a complete set of one hadron-states into the 
correlation function and we consider the single-particle K and states, this way we obtain 



< 


J n 


K{p) >< KK 


ao>< ao{p + q)\j''o 


> 




\p + qf-Ml 


(p2 - Ml) 



+ [ ^'[, r 2, ' (11) 

Jso Js'o [s - {p + q)^\ {s' - p^) 

where the hadronic spectral density includes the contributions of higher resonances and 
the hadronic continuum. The matrix element < KK \ oq > defines the coupling constant 

< K+{q)K-{p) I ao{p + q) >= QaoK+K- (12) 
and the current-particle matrix elements are given as 

<ao{p + q) \f°\0>^Xao , (13) 

< I I K{p) >= ifKP^ . (14) 
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After performing a similar double Borel transformation we obtain for the hadronic repre- 
sentation the result 



+ f ds f ds' p^''"\s,s')e 



so 



s/Ml-s'/Ml 



(15) 



The sum rule for the coupling constant gaoK+K- then follows by equating the expressions 
Ti(Mf,M|) obtained for the invariant function Ti{{p + q)'^,p'^) by theoretical (Eq. 8) and 
by physical (Eq. 15) considerations. In order to do this we have to identify the second 
term in Eq. 15 representing the continuum contribution with a part of the term calculated 
theoretically in QCD, thus affecting the subtraction of the continuum. The prescription that 
has been suggested for this purpose [12,13] is based on the observation that the distribution 
amplitudes (ppxiu) and (faxiu) are polynomials in {1 — u), therefore we can write 



N 



" A:=0 



(16) 



The continuum subtraction is affected in the leading twist 3 quark-antiquark term, since the 
contribution of the twist 3 quark-antiquark-gluon term in Eq. 8 is small, therefore finally 
we obtain the sum rule for the coupling constant gaoK+K- in the form 



9aoK+K- 



2A 



N 

Eh 

k=0 



Ml 



A' 



1=0 



JK Jo Juo- 



(p3K{oii,l - ai - a3,a3) [2 — 1]} (17) 



•ai as 



as 



where A — sq/M"^ with sq the smallest continuum threshold. 

In the numerical evaluation of the sum rule we use the twist 3 kaon distribution ampli- 
tudes given by [11] 



fpK{u) = 1 30 



3K 



+ 



JZKUJ^K 27 

* ^J^KfK 20 



1 + 6af ) 



'(2« - 1) 



(18) 



VaK{u) = 6UU I 1 + 



t- f3K 1 X 



20 //I 



(19) 



^3k{u) = 360q;iQ;2Q;3 



1 + ^(7a3 - 3) 



(20) 



where C!^{2u — 1) are the Gegenbauer polynomials. The overlap amphtude has been 
determined previously as A^q = (0.21 ± 0.05) GeV^ employing QCD sum rules method [8]. 
We also adopt the values at the renormahzation scale 1 GeV ms{lGeV) — 150 MeV, in 
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SU(3)/« limit f3KilGeV) = fs^ilGeV) = 0.0035 GeV^, u;3K{lGeV) = -2.88 and fx = 
0.160 GeV [11] with Mk = 0.4937 MeV. 

We then study the dependence of the sum rule for the coupling constant gaoK+K- on the 
continuum threshold sq ^-nd on the Borel parameters Mf and M| by considering independent 
variations of these parameters. We find that the sum rule is quite stable for the range of 
these parameters 1.00 < sq < 1.10 GeV^, 0.7 < < 1.4 GeV^ and 2 < M| < 6.0 GeV^. 
By varying the values of the parameters Sq, Mf, and M| in these regions we obtain the 
result for the coupling constant gaoK+K- as 4.4 < QaoK+K- < 5.6 GeV . The variation of 
the coupling constant as a function of the Borel parameters and M|, and the continuum 
threshold Sq is shown in Fig. 1. We note that the sign that we obtain for the couphng 
constant gaoK+K- is negative, that is gaoK+K- < 0. 

There has been several previous estimations of the coupling constant gaoK+K-- The 
Novosibirsk SND collaboration data of the radiative decay (f) n^rj^y [14] was analyzed 
in a phenomenological framework in which the contributions of p meson, chiral loop and 
Oo meson were considered, and the value gaoK+K- — (—1-5 ± 0.3) GeV was obtained for 
this coupling constant [15]. The coupling constant thus obtained results in constructive 
interference between the contribution of different amplitudes. From the analysis of their 
experimental data of ^ TT^r]'-/ decay, the KLOE collaboration estimated the coupling 
constant gaoK+K- as gaoK+K- = (2.3 ±0.7) GeV [16]. A new analysis of the KLOE data on 
(f) — > n^rj^ decay, on the other hand, gives the result gaoK+K- — {'^■^^-I'M) GeV [17]. In this 
analysis the phase S of the interference between — >• ao7 —>■ TT^rjj and p°7r° tt^t)^ 
amplitudes was obtained as 5 = 0, which is in accordance with the constructive interference 
of the different amplitudes observed in the phenomenological analysis of (I) ^ n^riy decay 
[14]. Our estimation of the coupling constant gaoK+K- using light cone QCD sum rules 
results in a value somewhat larger than the previous determinations based on the analysis of 
— > TT^T)^ data. However, because of the intrinsic uncertainties of the light cone QCD sum 
rule method, our result can be taken only to indicate that the scalar ao meson state may 
have somewhat large strong coupling. Finally wc note that the sign of the coupling constant 
gfoK+K- was obtained as gjQK+K- > in a previous light cone QCD sum rule determination 
of this couphng constant [13], and in a phenomenological analysis of — > 7r°7r°7 decay [18]. 
Thus our result about gaoK+K- that is gaoK+K- < 0, is consistent with the result obtained 
in the q'^q'^ model where g^K+K- = -QaoK+K- [2,5,19]. 

It has been argued that the ratio R = BR{(f) /07) / BR{(f) — > ao7) can provide insight 
into the structure of the scalar ao(980) and /o(980) mesons [4]. The KLOE collaboration 
obtained this ratio as [16] 

R^^^^:^^^6.1±0.6 . (21) 

This ratio can be written, assuming the intermediate KK loop mechanism for these decays, 
in the form [4] 

fl = |^iM-'^» (22) 

9aoK+K- ^aoV^) 

where the factors Fj^{R) and F^^{R) are related to the spatial extensions of /o(980) and 
ao(980) mesons, and for point-hke effective field theory calculations F^(i?) — 1. For a 
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spatially extended system with r.m.s. radius R > 0{Aqqjj) the high momentum region of 
the integration is suppressed [4], resulting in the form factor with the property F^{K) < 1, 
conversely F'^{R) — 1 means a spatially compact system. The angle 9 is the isospin mixing 
angle in the /o — ao system. If we use the result 6.2 < Qf^K+K- < 7.8 GeV obtained by 
a light cone QCD sum rule calculation [13], and our result 4.4 <| Qf^K^K- \< 5.6 GeV 
obtained by a similar light cone QCD sum rule method calculation we obtain 

Since in light cone QCD sum rule calculations isospin is assumed to be exact, which corre- 
sponds to ^ = 45°, we thus find FfQ(R)/Fao{R) ~ 1.8 which seems to imply that the spatial 
extensions of /o(980) and ao(980) mesons are not equal. Therefore, also given the relation 
that we obtain about the relative sign between the coupling constants gfyK+K- and gaoK+K- 
in accordance with g^g^ model, we may suggest that our result supports the view that the 
structure of ao(980) and /o(980) mesons is a combination of a KK molecule with a compact 
q^q'^ core with the spatial extension of /o(980) being more compact than that of ao(980) 
meson. 
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FIG. 1. The coupling constant gaoK+K- ^ ^ function of the Borel parameter for different 
values of the threshold parameter sq and the Borel parameter . The curves denote the limits 
of the stability region. 
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